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Abstract 

Open neutrino physics issues require precision studies, both theoretical and experimental ones, 
and towards this aim coherent neutral current neutrino-nucleus scattering events are expected to 
be observed soon. In this work, we explore n-nucleus processes from a nuclear theory point of view 
and obtain results with high confidence level based on accurate nuclear structure cross sections 
calculations. Besides cross sections, the present study includes simulated signals expected to be 
recorded by nuclear detectors, differential event rates as well as total number of events predicted 
to be measured. Our original cross sections calculations are focused on measurable rates for the 
Standard Model process, but we also perform calculations for various channels of the non-standard 
neutrino-nucleus reactions and come out with promising results within the current upper limits of 
the corresponding exotic parameters. We concentrate on the possibility of detecting (i) supernova 
neutrinos by using massive detectors like those of the GERDA and SuperCDMS dark matter 
experiments and (ii) laboratory neutrinos produced near the spallation neutron source facilities 
(at Oak Ridge National Lab) by the COHERENT experiment. Our nuclear calculations take 
advantage of the relevant experimental sensitivity and employ the severe bounds extracted for 
the exotic parameters entering the Lagrangians of various particle physics models and specifically 
those resulting from the charged lepton flavour violating —>■ e~ experiments (Mu2e and COMET 
experiments). 
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I. INTRODUCTION 


Coherent scattering of nentrinos on complex nnclei was proposed long ago [i|, as a 
prominent probe to study neutral-current (NC) i/-nucleus processes, but up to now no events 
have been experimentally measured. Neutrino detection, constitutes an excellent probe to 
search for a plethora of conventional neutrino physics applications and new-physics open 
issues js-b]. In principle, low-energy astrophysical and laboratory neutrino searches provide 
crucial information towards understanding the underling physics of the fundamental elec- 


B- 


troweak interactions within and beyond the SM l6ll. Well-known neutrino sources include fi) 
supernova neutrinos (with energies up to 60-100 MeV) and (ii) laboratory neutrinos (with 
energies up to 52.8 MeV) emerging from stopped-pion and muon decays at muon factories 


(Fermilab, PSI, JPARC, etc.) and at the spal 


ation neutron source (SNS) at Oak Ridge Na¬ 


tional Lab Recently, it became feasible jsl to detect neutrinos by exploiting the neutral 
current interactions and measuring the nuclear recoil signals through the use of very low 
threshold-energy detectors [ol, ll0|. To this purpose, great experimental effort has been put 
and new experiments have been proposed to be performed at facilities with stopped-pion 
neutrino beams, based on promising nuclear detectors like those of the COHERENT experi¬ 
ment III, ll2j and others Il3l at the SNS, or alternative setups at the Booster Neutrino Beam 
(BNB) at Fermilab Q 0. The ..clear .-detectore adopted by the relevant experiments 
include liquid noble gases, such as ^°Ne, "^^Ar, ^^^Xe as well as, ^®Ge and CsI[Na] detection 


materials 


16|- 


On the theoretical side, the t/-signals of low-energy neutrinos, expected to be recorded 
in sensitive nuclear detectors [l7, 18 1, could be simulated through nuclear calculations of 
i/-nucleus scattering cross sections. Such results may provide useful information relevant for 
the evolution of distant stars, the core collapse supernovae, explosive nucleosynthesis and 
other phenomena 2, 201. In fact, coherent neutral current i/-nucleus scatterin g event s are 


expected to be observed by using the high intensity stopped-pion neutrino beams 


21 


221 and 


nuclear targets for which recoil energies are of the order of a few to tens of keV, and therefore 


appropriate for detection of WIMPs 
detectors are e.g. the SuperCDMS 


24| |. candic 




27|, GERDA 


ates of cold dark-matter 


25 


261. Such 


and other multi-purpose detectors 


29N31|. For low-energies, the dominant vector components of NG interactions 


coherent contribution of all nucleons (actually all neutrons) in the target nucleus 


eac to a 


32 


34|. 
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It is worth mentioning that, after the discovery 35l-l39| of neutrino oscillations in propa¬ 


gation, the challenge of neutral and chargee 
tigated by extremely sensitive experiments 


l^ton flavour violation (LFV) is further inves- 


40l-l44| searching for physics beyond the current 


Standard Model (SM) To this end, neutrino-nucleus coherent scattering experiments 


may pro 

0 , 3-48 


3e new physics beyond the SM involved in exotic neutrino-nucleus interactions 


48|, an undoubtable signature of non-standard physics. Therefore, new data and in¬ 
sights will be provided to the physics of flavour changing neutral current (FCNC) processes. 


in the leptonic sector, in non-standard neutrino oscillation effects 491 451|, in neutrino tran¬ 


sition magnetic moments 


such experimental sensitivity may a. 
open issues in nuclear astrophysics 
In recent works 


52| . in sterile neutrino search 53|] and others 5^. Furthermore, 
so inspire advantageous probes to shed light on various 


55 


56| 


48j, neutral-current (NC) non-standard interactions (NSI) involving 


(anti)neutrino scattering processes on leptons, nucleons and nuclei have been investigated. 
The reactions of this type that take place in nuclei are represented by 


Z) —)■ Vp) -|- (A, Z) , 


( 1 ) 


(a, (3 = e, with a ^ (3). It has been suggested [57| that, theoretically the latter processes, 
can be studied with the same nuclear methods as the exotic cLFV process of /i“ e~ 


conversion in nuclei 


58 


sol. The eon-esponding Lageapgians may be denved w.thin the 


context of various extensions of the SM [a, l6ll], like 


seesaw model 


62 


63| . left-right symmetric models 6^, gluonic operator model 65j, etc. 


he four fermion contact interaction. 


It is well-known that neutrino NSI may have rather signiheant impact in many areas of 
modern physics research and thus, motivate a great number of similar studies 66|]. Especially 
in astrophysical applications, constraints coming out of some supernova explosion scenarios 


67 


69|, may be affected and eventually lead to the necessity of further investigation of 


NSI in both LFV and cLFV processes that may occur in solar and supernova environment 


70 


72| . Such open issues motivated our present work too. 


One of our main purposes in this paper, which is an extension of our previous study 48| . 
is to comprehensively study the above issues by performing nuclear structure calculations 
for a set of experimentally interesting nuclei. We estimate reliably the nuclear matrix ele¬ 
ments describing both interaction channels, the exotic and the Standard Model ones, but 
we mainly focus on the SM component of the neutrino-nucleus processes, i.e. we consider 
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a = /3 in the reactions of Eq. ©■ Exotic neutrino-nucleus events are also computed. By ex¬ 
ploiting our accurate original cross sections, we obtain simulated //-signals and flux averaged 
cross sections which are experimentally interesting quantities for both Supernova and SNS 
neutrinos. The total number of events expected to be recorded over the energy-threshold 
for the studied nuclear targets are also presented for both cases. 

We stress that, we have devoted special effort to obtain results of high accuracy by 
constructing the nuclear ground state within the context of the quasi-particle random phase 
approximation (QRPA), i.e. by solving iteratively the BCS equations for realistic pairing 


interactions (the Bonn C-D pot ential) 73l. l74j . and achieving high reproducibility of the 
available experimental data In addition, we made comparisons with the results of 


other methods evaluating the nuclear form factors that enter the coherent rate 76|, l77| as 


the one which employs fractional occupation probabilities (FOP) of the states (on the basis 
of analytic expressions) 78|, and other well-known methods 79 |. 


II. DESCRIPTION OF THE FORMALISM 

In this section we present briefly the necessary formalism for describing all channels of 
the NSI processes of the reactions ([I]), derived by starting from the corresponding nuclear- 
level Feynman diagrams. In Fig. [T] the exchange of a Z-boson between a lepton and a 



FIG. 1. Nuclear level Feynman diagrams for: (a) SM Z-exchange neutral current //-nucleus reac¬ 
tions, (6) non-standard Z-exchange //-nucleus reactions, and (c) Z-exchange and photon-exchange 
—>■ e~ in the presence of a nucleus (muon-to-electron conversion). The non-standard (cLFV or 
LFV) physics enters in the complicated vertex denoted by the bullet • 48l |. 
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nucleon is represented, for the SM i/-nucleus scattering, Fig. [I](a), and for the exotic z/- 
nucleus scattering. Fig. Ht^b). As already mentioned in the Introduction, the non-standard 


z/-nucleus processes and the exotic cLFV /i“ e~ conversion in nuclei [2, H, 7^, 
can be predicted within the context of the same new-physics models j^, For this reason, 
in Fig. dl^c) we also show the exchange of a Z-boson or a virtual 7 -photon leading to the 


nuclear fi e conversion, [^,1^. Thus, the leptonic vertex in the cases of Fig. [I]^b),(c) 


is a complicated one. A general effective Lagrangian that involves SM interactions (£sm) 
and NSI (/Insi) with a non-standard flavour preserving (FP) term, a non-universal (NU) 
term and a flavour changing (FC) term reads 




SM 


-^NSI — + C. 


-NU 


£ 


FC • 


( 2 ) 


Each of the components £sm and £nsi, the individual terms and Lpc as well as the 
nuclear matrix elements that arise from each part, are discussed below. 

A. Coherent cross sections of non-standard u-nucleus reactions 


The quark-level Lagrangian for neutral current non-standard neutrino interactions £nsI) 
at the four fermion approximation, (energies -C Mz) is parametrized as j^, 67| 


£nsi = -2V2GF [fYPf] , 


(3) 


f=u,d 

a,p=e,ii,T 


where / denotes a hrst generation SM quark, are three light neutrinos with Majorana 
masses and P = {L,R} are the chiral projectors. In the latter Lagrangian (j3]), two classes 
of non-standard terms are considered (i) flavour preserving non-SM terms that are pro¬ 
portional to (known as non-universal, NU interactions) and (ii) flavour-changing (FC) 
terms proportional to a ^ (5. These couplings are dehned with respect to the strength 
of the Fermi coupling constant Gi 


47 


67l | . In the present work we examine spin-zero nuclei 
thus, the polar-vector couplings dehned as -f- are mainly of interest. For the 

axial-vector couplings it holds 


* 


Following Ref. 72|, the nuclear physics aspects of the neutrino-matter NSI can be ex¬ 
plored by transforming the quark-level Lagrangian ([3]), eventually to the nuclear level where 
the hadronic current is written in terms of NC nucleon form factors that are functions of 
the four momentum transfer. Generally, for inelastic u-nucleus scattering, the magnitude of 
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the three momentum transfer, g = |g|, is a function of the scattering angle of the outgoing 
neutrino 6 (in laboratory frame), the initial, Ei, and hnal, Ef, nuclear energies, as well as 
the excitation energy of the target nucleus, uj, and takes the form = u‘^ + 2EiEf (1 — cos 6) 
[^, 76|. Our analysis in the present paper, concentrates on the dominant coherent (elastic) 
channel where only gs —?• gs transitions occur {uj = 0, = Ef) and the momentum transfer 

in terms of the incoming neutrino energy, E,y, becomes g^ = 2E‘^{1 — cos6) or equivalently 
g = 2Ei, sin(6'/2). 


The NSI coherent differential cross section o: 


neutrinos scattering off a spin-zero nucleus. 


with respect to the scattering angle 9 reads [48 1 




( 4 ) 


a = e,g.,T, denotes the flavour of incident neutrinos and |gs) represents the nuclear ground 


state (for even-even nuclei assumed here, \gs) = | J"") = 


that enters the cross section of Eq. (jlj), is written as 481 


0^)). The nuclear matrix element. 


\M 


NSI |2 _ 
V.Va 


^\{9s\\G-y^im9s)\ = 

ZE,{q-) + iVF^(g^)]^ 

E ZFz{q^) + {efp + 2e^^) iVF^(g^)] ^ , 


( 5 ) 


(/5 = e,fi,T) where Ez(n) denote the nuclear (electromagnetic) form factors for protons 
(neutrons). We stress that, in the adopted NSI model, the coherent NC i/-nucleus cross 
section is not flavour blind as in the SM case. Obviously, by incorporating the nuclear 
structure details, in Eqs. (jl]) and ([S]), the cross sections become more realistic and accurate 
js]. The structure of the Lagrangian (|5]), implies that in the r.h.s. of Eq. ([S]) the hrst term 
is the NU matrix element, and the summation is the EC matrix element, 

hence we write 




V,u^\ 




( 6 ) 


From experimental physics perspectives, it is rather crucial to express the differential cross 
section with respect to the recoil energy of the nuclear target, T^. In recent years it became 
feasible for terrestrial neutrino detectors to detect neutrino events by measuring nuclear 
recoil [^, 16|. Therefore, it is important to compute also the differential cross sections 
da/dTjy. In the coherent process, the nucleus recoils (intrinsically it remains unchanged) 
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with energy which, in the approximation ^ takes the maximum value = 

2El/[M + 2Ejj), with M denoting the nuclear mass |32|, l33|. Then, to a good approximation, 
the square of the three momentum transfer, is equal to = 2MT;v, and the coherent NSI 
differential cross section with respect to T/v can be cast in the form 

_ GjpM f ^ MTm 


1 - 


(ITn tt 

We note that, compared to previous studies 




|(^s||G'y®i(g)||^s)| 


(7) 


55 


66l |. we have also taken into consideration 


the interaction z/ — m quark [see Eq. ([5])], in addition to the momentum dependence of the 


nuclear form factors 


^]. Both Eqs. 


) and ([7]) are useful for studying the nuclear physics 

of NSI of neutrinos with matter. 

Furthermore, by performing numerical integrations to Eq. (jl]) over the scattering angle 6 
or to Eq. ([7|) over the recoil energy Tjv, one can obtain integrated (total) coherent NSI cross 
sections, ctnsi,i/„- Following Eq. ([S]), the individual cross sections crNu.^a (Tpc,ua be 
evaluated accordingly. 


B. Standard Model coherent i/-nucleus cross sections 


The effective (quark-level) SM z/-nucleus interaction Lagrangian, £sm at low and inter¬ 
mediate neutrino energies, is written as 


= -2V2Gf ^ Qp [lyalpLuo] [fYPf] , 


( 8 ) 


/= u,d 

a=e,ii,T 


where g'f^ = ^ § sin^ 9w and gp = sin^ 9w are the left- and right-handed couplings of the 

M-quark to the Z-boson and S'l = — | | sin^ 9w and gp = 4 sin ^ 9w are the corresponding 

couplings of the d-quark {9^ is the Weinberg mixing angle) [76|. 

For coherent z/-nucleus scattering, the SM angle-differential cross section reads 


= ^EI{1 + cos 9) {gs\\Mo{q)\\gs) 
a cos 9 271 


(9) 


The operator j\4o in the nuclear matrix element of the latter equation is the Coulomb 
operator which is equal to the product of the zero-order spherical Bessel function times the 
zero-order spherical harmonic |73|, l76|. This matrix element can be cast in the form 711 

2 


\M 


SM |2 _ 
V,Uc = 


{9s||M„(«)||(,s> = [fvZFzif) + g^NF,t(q'‘)Y , 


( 10 ) 


7 





















where, the polar-vector couplings of protons gy and neutrons Qy with the Z boson (see Fig. 
[IKa)), are written = 2{gl +g'^) + {gi +gj^) = l-2sm^ew and g^ = {gl +g^)+ 2{g'l + 
gjf) = — |, respectively. As can be easily seen, the vector contribution of all protons is very 
small {gy ~ 0.04), hence the coherence in Eq. (ITOil essentially refers to all neutrons only 
of the studied nucleus. After some straightforward elaboration the differential cross section 
with respect to the nuclear recoil energy, T/v, takes the form 


daf 


SMji/cj 


dTi 


N 


GlM 

71 


MT, 


N 


2El 


(^s||A4o(g)||^s) 


( 11 ) 


The Lagrangian £tot of Eq. ([2]), contains the flavour preserving (FP) part, equal to 
£fp = which can be evaluated through the Coulomb matrix element 




SM 


= 


+ M 


NU 

V,Ua 


( 12 ) 


Subsequently, the total coherent cross section may be computed on the basis of the matrix 
element 

. r» __ o _ . O 

(13) 


|X“. 


r = |viyt 


+ |nE, 


In a previous work we evaluated original differential cross sections dax^u^/d cos 9 and 
dax^^^/dT]\f, as well as individual angle-integrated cross sections of the form ax,i,^{E^), with 
a = and A = SM,NU,FP,FC (FC stands for the six flavour changing processes 

In this work, we perform Standard model cross sections calculations (for convenience 
from now on we drop the index A = SM and always consider for a set of nuclei 

throughout the periodic table up to ^°®Pb. We adopt various nuclear models (see section 
HID to compute the nuclear form factors. Then, for a great part of the cross section results 
(except differential cross sections) we evaluate folded cross sections, and event rates. 


III. EVALUATION OF THE NUCLEAR FORM FACTORS 

A. Nuclear Structure calculations 

At hrst, we study the nuclear structure details of the matrix elements entering Eq. flTOl) . 
such results reflect the dependence of the coherent cross section on the incident-neutrino 
energy Ey and the scattering angle 9 (or the recoil energy Tjv). We mention that for the 
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even-even nuclei this study involves realistic QRPA calculations for the differential cross 
sections dau^/dcos9 and da^^/dT^, performed after constructing the nuclear ground state 
\gs) by solving iteratively the Bardeen Cooper Schriffer (BCS) equations. The solution of 
these equations provides the probability amplitudes and of the j-th single nucleon 
level to be occupied or unoccupied, respectively. Moreover, the latter equations provide the 
single quasi-particle energies, based on the single particle energies of the nuclear held (a 
Coulomb corrected Woods-Saxon potential in our case) as well as the pairing part of the 
residual two-body interaction (Bonn C-D potential in our case). Then, the nuclear form 
factors for protons (neutrons) are obtained as 711 


N„ 




(14) 


with [j] = \/2j -|- 1, Nn = Z (or N). For each nuclear system studied, the chosen active 
model space, the harmonic oscillator (h.o.) parameter b and the values of the two parameters 
9p!fr proton (neutron) pairs that renormalise the monopole (pairing) residual interaction 


(obtained from the Bonn C-D two-body potential describing the strong two-nucleon forces), 
are presented in Table [H The a djust ment of is achieved through the reproducibility of 
the pairing gaps Ap(„) (see e.g. j20 | ). 


TABLE I. The values of proton and neutron pairs that renormalise the residual interaction 
and reproduce the respective empirical pairing gaps Ap and A^. The active model space and the 


harmonic oscillator parameter, for each isotope, are also presented. 

Nucleus 

model-space b Ap g^^^^ 

12c 

8 (no core) 1.522 4.68536 4.84431 1.12890 1.19648 

160 

8 (no core) 1.675 3.36181 3.49040 1.06981 1.13636 

20Ne 

10 (no core) 1.727 3.81516 3.83313 1.15397 1.27600 

28 si 

10 (no core) 1.809 3.03777 3.14277 1.15568 1.23135 

40 Ar 

15 (no core) 1.902 1.75518 1.76002 0.94388 1.01348 

48Ti 

15 (no core) 1.952 1.91109 1.55733 1.05640 0.99890 

™Ge 

15 (no core) 2.086 1.52130 1.56935 0.95166 1.17774 

i32Xe 

15 (core ^°Ca) 2.262 1.19766 1.20823 0.98207 1.13370 


9 











B. Other methods for obtaining the nuclear form factors 


The nuclear form factor, which is the Fourier transform of the nuclear charge density 
distribution Pp{r), is dehned as 

47r 


Fz{q^) = 


Pp{r)jo{qr)r dr, 


(15) 


with jo being the zero order spherical Bessel function. Due to the signihcance of the nuclear 
form factors in our calculations and for the benefit of the reader we devote a separate 
discussion to summarise some useful possibilities of obtaining these observables. 

i) Use of available experimental data 

For many nuclei and especially for odd-A isotopes, the proton nuclear form factors Fz{q^), 
are computed by means of a model independent analysis (using a Fourier-Bessel expansion 
model or others) of the electron scattering data for the proton charge density Pp{r) [75| 
wherever, possible. The absence of similar data for neutron densities, restricts us to assume 
that F^iq^) = Fz{q^)- In this work, we consider this method only for the case of the very 
heavy doubly closed nucleus. 

ii) Fractional occupation probabilities in a simple Shell-Model 


In Ref. 


78| the form factor Fz{q^), for h.o. wavefunctions has been written as 701 






4(«6,Z)= 


( 16 ) 


A=0 


The radial nuclear charge density distribution Pp{r), entering the dehnition of Eq. (1T5|1 . is 
written in the following compact form 701 


-| ^max 

PpM = n (I z) , n (X, z) = ^ 


(17) 


A=0 


where y = r /b, with b denoting the h.o. size parameter. A"max = (2u -|- f')max stands for the 
number of quanta of the highest occupied proton (neutron) level. The coefficients f\ are 
expressed as 

_ MU (On -r 

( 18 ) 


^ + l)n! ^ 


K 6 . + ’ 

where F(^is the Gamma function. For the coefficients 9\, and further information 


see Ref. 


7Q|. 
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FIG. 2. The charge density distribution (left) and the form factor as a function of the momentum 
transfer {right), for the cases of ^*^Ar and ^®Ti nuclei. The introduction of fractional occupation 
probabilities (FOP) of the states provides higher reproducibility of the experimental data, compared 
to the simple Shell-Model and that of Eq. (I2ip . The BCS nuclear neutron form factor Fj^{q^) is 
also presented and compared. 


Up to this point, that the proton occupation probabilities entering Eq. ffTHll Eq. ffTHll have 
been considered equal to unity for the states below the Fermi surface and zero for those above 


the Fermi surface. In Ref. 


78| , the authors introduced depletion and occupation numbers, to 
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parametrise the partially occupied levels of the states. These parameters satisfy the relation 


+ 1 ) — 


(19) 


{n£)j 

all 


Within this context, the ’’active” surface nucleons (above or below the Fermi level) have 
non-zero occupation probability anej 7 ^ 0, smaller than unity, while the ’’core” levels have 
occupation probability anij = 1. In this paper we extend the work of Ref. 78| where three 
parameters ai, 02 , are used to describe the partial occupation probabilities of the surface 
orbits. We improve the formalism by introducing more parameters, increasing this way 


the number of ’’active” nucleons in the stuc 


ied nuclear system and come out with higher 


reproducibility of the experimental data 7^. To this aim, we introduce four parameters 


Q!j, i = 1,2,3,4 in Eq. ([19]). Then, the assumed ’’active” single-particle levels are hve and 
Eq. (16) of Ref. Q becomes 

^{x,Z,ai) =n(x,Z2) 


ai 


+ n(x, Zc) 
+ n(x,z') 
+ n(x,z"; 


Zi — Z 2 

Z' -z 


n(x,^i 


02 




OL 2 


Z' -z, 
z-z. 


+ 


Zc — Z\ 

«3 


Z' - Zc Z' - z. 


04 


A 


as 

Z' -Zc 

a4 

Z" - Z' 

+ n(x,z'") 


( 20 ) 


A 


Z'" - Z'' 


Z” - Z' Z'" - Z" 

with A = ai -|- a 2 — as — 04 . By substituting the polynomial n(x, Z) of Eq. flT7|l with 
that of the latter expression and using the experimental data 75| , we we £t the parameters 
ttj (and similarly for the form factor of Eq. ffT 6 |) ). As an example, for the ^°Ar isotope 
we have, Z 2 = 10, Z^ = 12, Z = Z,. = 18, Z' = 20, Z" = 22, Z'" = 30. The resulting 
fractional occupation probabilities that £t the experimental charge density distribution are 
ai = 0.85, a 2 = 1.25, = 0.85, a 4 = 0.75. Similarly for the "^®Ti nucleus, we have 

Z 2 = 18, Zi = 20, Z = Zc = 22, Z' = 30, Z" = 34, Z'" = 40 and the fitting parameters 
are ai = 1.0, 02 = 1.5, = 0.35, a 4 = 0.1. In Fig. [2] the prediction of the method is 

compared with that of the simple shell-model and the experimental data. We note that in 
the momentum transfer range of our interest (i.e. q < 2 fm“^) the form factor has excellent 
behaviour. We however mention that even though the FOP method presents very high 
reproducibility of the experimental data, it is not always applicable, e.g. for deformed nuclei 
(where BCS appears to be still successful). 


12 




























iii) Use of effective expressions for the nuclear form factors 
We finally discuss one of the most accurate effective methods for calculating the nuclear 
form factor by Ref. 791 

r 1 J 

( 21 ) 


F{q ) = - - -exp 


qRo 




where ji{x) is the known hrst-order Spherical-Bessel function and Rq = — 5s^, with R 

and s being the radius and surface thickness parameters of the nucleus respectively. The 
radius parameter is usually given from the semi-empirical form R = 1.2^4^/^ fm while s is of 
the order of 0.5 fm (see Ref. [^). 

It is worth noting that, by inserting the form factors Fz{n) obtained as described above 
in Eq. (fT0|) . the resulting cross sections have a rather high conhdence level. In the next part 
of the paper the results show that the momentum dependence of the nuclear form factors 
becomes crucial, especially for intermediate and high energies. In some cases, differences 
of even an order of magnitude may occur as compared to the calculations neglecting the 
momentum dependence of the nuclear form factors. 


IV. RESULTS AND DISCUSSION 

A. Integrated coherent n-nuclens cross sections 

The next phase of our calculational procedure is related to the total coherent z/-nucleus 
cross sections, obtained through numerical integration of Eq. ([9]) over angles [or Eq. CD 
over T^r] as 

^uAEu) = [ -p^Ae.E,) dcose. (22) 

J a cos 

The results for the Standard Model cross sections, for a set of different promising targets 
throughout the periodic table, are presented in Fig. |3l As can be seen, the present nuclear 
structure calculations indicate that between light and heavy nuclear systems, the cross 
sections may differ by even two orders of magnitude (or more) as a consequence of the 
dependence on the nuclear parameters (i.e. mass, form factors, etc.). We also see that for 
heavier nuclei the cross sections flatten more quickly (at lower neutrino energies) compared 
to that of lighter nuclear isotopes. The latter conclusion originates mainly from the fact 
that, for heavy nuclei the suppression of the cross sections due to the nuclear form factors 
becomes more signihcant. Thus, for heavy material the nuclear effects become important 
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even at low energies. Such original cross section results are helpful for the simulations of 
the Standard and non-standard Model signals of z/-detection experiments (see below). 



FIG. 3. Total coherent cross sections in units 10 ^®cm^ for a set of nuclei as a function of 

the incoming neutrino energy Ei,, for the SM neutrino processes + {A, Z) —)• Vai^a) + (d., Z). 


B. Supernova neutrino simulations 


As discussed previously, our present calculations may also be useful for ongoing and 
future neutrino experiments related to supernova (SN) neutrino detection, since as it is 
known, the neutrinos emitted in SN explosions transfer the maximum part of the total the 
energy released. Then, the total neutrino flux, arriving at a terrestrial detector as 


a function of the SN neutrino energy Ey, the number of emitted 


distance d from the source (here we consider d = lOkpc), reads 23 


anti)neutrinos at a 


34| 


Ny 




(23) 


(a = e,/r,r) where denotes the energy distribution of the (anti)neutrino flavour a. 

The emitted SN-neutrino energy spectra rj^^[Ey) may be parametrised by Maxwell- 
Boltzmann distributions that depend only on the temperature Ty^ of the (anti)neutrino 
flavour Va or z/„ (the chemical potential is ignored) we have 


'lua 2T3 

U<x 


( 24 ) 
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= 3.5MeV, = S.OMeV, = S.OMeV, x = /i,r 32|). For each flavour, the total 


number of emitted neutrinos Ny^ is obtained from the mean neutrino energy 


48 | 




(25) 


and the total energy released from a SN explosion, U = 3 x 10^'^erg 


17 |. 


C. Laboratory neutrino simulations 

The spallation neutron source (SNS) at Oak Ridge National Lab produces neutrons 
by firing a pulsed proton beam at a liquid mercury target 54l. The main aim of the 


COHERENT proposal 


11 


, ll2| (or of other similar conceptual 


13 


14| 1 concerns with possible 


detection of neutrino-nucleus coherent scattering events at the SNS. Our simulations here, 
are mainly motivated by previous studies js, 15, 16, 53| and the hope to provide our accurate 
nuclear structure calculations. 

In stopped pion-muon sources, neutrinos are produced by the pion decay chain. Pion 
decay at rest 7r+ —)■ (r = 26 ns) produces monochromatic muon neutrinos at 29.9 

MeV, followed by electron neutrinos z/g and muon antineutrinos that are produced by the 
muon-decay /r’*' —)■ (r = 2.2 /rs) [21 


22l |. For pulsed beams in time-scales narrower 
than /is, z/g’s and h^’s will be delayed with the beam while z/^’s will be prompt with the beam 
Sj. The emitted Ve and i)^ neutrino spectra are described by the high precision normalized 
distributions, known as the Michel spectrum 


=9QElM;^ {M, - 2Ey) , 
=lQElM-^ {3M, - AEy) , 


(26) 


(M^ = 105.6MeV is the muon rest mass). The maximum neutrino energy in the latter 
distributions is = M^/2 = 52.8MeV (see e.g. 9]). 

The spallation neutron source (SNS) at Oak Ridge National Lab is currently the most 
powerful facility to detect for a first time neutrino-nucleus coherent scattering events, since 
it provides exceptionally intense fluxes = 2.5 x 10^ z/s“^cm^ at 20 m and <F,y„ = 6.3 x 


10® z/s ^cm^ at 40 m from the source 


21 


22| . The simulated laboratory neutrino signals 


a 


Sign 

u,lab. 


coming out of our calculations for the adopted nuclear targets are discussed below. 
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D. Simulated neutrino signals 


By weighting the integrated cross section av^{Ey) with the neutrino distributions of Eq. 
), for SN neutrinos, or Eq. fl2^ . for laboratory neutrinos, the total signal produced on a 
terrestrial detector is described by [74 1 


= ^ = SN,iat. (27) 

The resulting signals, cr**|"'(Ej,), obtained by inserting in Eq. fl27|) the cross sections of 
Fig. El are plotted in Fig. m 



FIG. 4. The signal cross sections that represent the expected signal to be recorded on a terrestrial 
nuclear j/-detector, (left) for supernova neutrinos = SN), evaluated with Maxwell-Boltzmann 
distributions at d = lOkpc, and (right) for SNS neutrinos = lab.), at 20 m from the source. For 
the case of SNS neutrinos the figure takes into account only the delayed beam, evaluated with the 
generic flux of ~ 10^ i 2 s“^cm“^. Different nuclear detectors have been studied. 


In our previous work 48| it was shown that the simulated cross sections reflect the char¬ 
acteristics of the incident neutrino spectrum of the specihc neutrino flavour a and therefore, 
such a simulated signal is characterised by its own position of the maximum peak and width 
of the distribution We, however, remind that within the framework of the SM, co¬ 

herent neutrino scattering is a flavour blind and a particle-antiparticle blind process. For 
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this particular case our results are shown in Fig. 0] for supernova and laboratory (SNS) 
neutrinos. 

In neutrino simulations, another useful quantity is the flux averaged cross section 
which in our notation is written as 

V'a 

The results for (cXy)^, obtained by using the angle-integrated cross sections of Fig. 0] are 


TABLE II. Flux averaged cross sections in units 10“^°cm^ for the adopted supernova {d = 

lOkpc) and laboratory (delayed flux only) neutrino spectra. For the case of SNS neutrinos, we 
adopt the generic flux, i.e. ~ z^s“^cm^ at 20 m for all nuclear targets. 


Nucleus 


(O'^)SN 


12q 160 20Ne 28Si 40^1. 48^1 56pg 76q^ 208pp 


1.46 2.51 3.91 7.52 18.59 25.43 33.29 70.63 207.56 514.93 
3.07 5.33 8.13 15.52 37.91 51.50 67.02 139.83 395.59 949.50 


listed in Table [TTl for both neutrino sources. 


E. Differential and total event rates 


From experimental physics 


perspectives, predictions for the differential event rate, 


of a //-detector are crucial 23|. The usual expression for computing the yield in events is 


based on the neutrino flux, and is defined as [34 1 


dN 


u.(r„) = ^ = Ky ]<!>,, 14, dE, 


dau 


Ua 


dcos9 


S{Tn- 


2M 


d cos 9 , 


(29) 


where K = Ntargdtot. accounts for the total number of nuclei (atoms) in the detector material 
Ntarg. times the total time of exposure ttot- Using the latter equation, one concludes that, the 
lower the energy recoil, the larger the potentially detected number of events (see Fig. [5] and 
Fig. |6]). In principle, in order to maximize the potential detection of a rare event process like 
the //-nucleus scattering, detector materials with very low energy-recoil threshold and low- 
background are required. In the last stage of our study we make predictions for the total 
number of coherent scattering events, the most important quantity, both from theoretical 
and experimental perspectives. To this purpose, we evaluate the number of expected counts. 
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FIG. 5. Yield in events (left) and total number of events over nuclear recoil threshold {right), 

for supernova neutrinos at d = 10 kpc. Here, 1 ton of perfectly efficient ^^Ne and ^^Ar detectors 
have been considered and also possible neutrino oscillation in propagation effects are neglected. 
For heavier nuclear targets the differences become rather significant. In this figure, F{q^) stands 
for Eq. (j2ip and FOP for the method of fractional occupation probabilities of the states. For more 
details see the text. 


for the studied detector materials, by performing numerical integration of Eq. (1291) over the 
nuclear recoil threshold (see Table [ITT|) . 


As has been discussed previously 


23 


24 . SN neutrino detection might become possible 


by the massive dark matter detectors [29| which have very good energy resolution and 
low threshold capabilities 34- These experiments are designed (or planned) to search for 


WIMPs 


m 


26| and/or other rare events such as the neutrinoless double beta decay. The 
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FIG. 6. Same as Fig. [5] but for and 


latter, use heavy nuclei as nuclear detectors, e.g. Ge (GERDA [28| and SuperGDMS [27 1 
experiments). In addition we report that SN neutrino events can be potentially detected by 
experiments using noble gases like Ne (CLEAN detector 2^), Ar (WARP programme 30l |) 
and Xe (XENON 100 Collaboration 31|h 


As mentioned in section HTTl in order to test our nuclear calculations we have also employed 
other nuclear methods. To this purpose, we have compared our original results evaluated 
with the BCS method with those obtained as discussed in subsection IIIIBI and concluded 
that for the case of the coherent channel all available nuclear methods are in good agreement, 
but their results differ signihcantly from those obtained assuming Fz{q^) = F^iq^) = 1 (see 
Fig. |5] and Fig. [ 6 ]). We stress however, that, since the cross section is mostly sensitive 
to the neutron distribution of the target nucleus, the most accurate method (at low and 
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TABLE III. Total number of events per ton of the target materials for a supernova at a distance 


of 10 kpc. We assume various energy t 
excellent agreement with those of Refs. 


ires 

12 

3 , 

34 


ds 5, 10, 25 or 50 keV. Our present results are in 


Nucleus 

Tn 

Tat > 5 keV T^r > 10 keV T^r > 25 keV T^r > 50 keV 

12c 

2.52 

2.25 

2.05 

1.60 

1.14 

160 

3.29 

2.84 

2.51 

1.83 

1.19 

20Ne 

4.03 

3.35 

2.87 

1.96 

1.16 

40 Ar 

9.46 

6.63 

5.01 

2.53 

1.00 

48Ti 

10.73 

7.04 

5.06 

2.27 

0.76 

56pe 

12.00 

7.36 

5.04 

2.01 

0.57 

76Ge 

18.58 

9.61 

5.82 

1.70 

0.30 

i32Xe 

30.68 

9.84 

4.16 

0.46 

0.01 

208pp 

46.93 

7.86 

1.95 

0.03 

< 10-3 


intermediate energies) is the BCS method which provides realistic proton as well as neutron 
form factors. All other methods employed here consider only the proton distribution and 
assume Fz{q^) = which especially for heavy nuclei, is a rather crude approximation. 

We remark, however, that the aforementioned nuclear methods offer reliable results on the 
differential and total event rates for low energies (see Fig. [5] and Fig. |6]), but in order 
to correctly estimate the neutron form factor, methods like the BSC are probably more 
appropriate. 

Our present nuclear structure calculations for laboratory (SNS) neutrinos (see Fig. 
ED, are in good agreement with previous results j^. They imply that a comparably large 
number of coherent neutrino scattering events is expected to be measured by using LNe, LAr, 
LXe, Ge and CsI[Na] materials adopted by the COHERENT Collaboration [n, 1^. The 
predictions of the BCS method for these nuclei are illustrated in Fig. [7] and compared with 
those of other promising nuclear targets. Because the neutrino flux produced at the SNS is 
very high, (of the order of ~ 10 ^i/s“^cm“^ per flavour at 20 m from the source 21 |). 
even kg-scale experiments expect to measure neutrino-nucleus coherent scattering events at 
signihcantly higher rates than those of supernova neutrinos. 

It is worth noting that, the choice of the target nucleus plays also a crucial role, since 
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a light nuclear target may yield almost constant number of events throughout the energy 
range, but small number of counts. On the other hand, a heavy nuclear target provides more 
counts, but yields low-energy recoil making the detection more difficult. This leads to the 
conclusion that the most appropriate choice for a nuclear detector might be a combination 
of light and heavy nuclear isotopes, like the scintillation detectors discussed in Ref. 341. 





FIG. 7. Total number of expected events over nuclear recoil threshold for 1 ton of various nuclear 
targets at 20 m from the source {^ua ~ 10^ i/s“^cm“^). The upper (lower) panel assumes the 
delayed (prompt) flux of laboratory stopped-pion neutrino sources. This figure assumes a perfectly 
efficient detector and negligible neutrino oscillation effects. 


21 
























F. Non-standard neutrino interactions at the COHERENT detector 


The multi-target approach of the COHERENT experiment 


11 


detection can also explore non-standard physics issues such as NSI 
moment 


47 


1211 aiming on neutrino 
El , neutrino magnetic 


52[ | and sterile neutrino j53|. In this subsection we find it interesting to evaluate the 


non-standard neutrino-nucleus events that could be potentially detected by this experiment 


B 


and 


in each of the proposed nuclear targets. The high intensity SNS neutrino beams 
the two promising //-detectors, liquid ^°Ne (391 kg) and liquid ^‘^Ar (456 kg) j^, hrstly 
proposed by the CLEAR [l3j and CLEAN [29(] designs (located at distance 20 m from the 
source), constitute excellent probes to search for the exotic //-reactions. Other possibilities 
111 . Il2l | include medium and heavy weight targets like ^®Ge (100 kg) inspired by the dark 
matter SuperCDMS 27] detector (located at 20 m) and ^^^Xe (100 kg located at 40 m). 

In Fig. [H]and Fig. [Hlthe resulting number of exotic events are illustrated and compared 
with the SM predictions. We note, however, that especially for the case of the flavour 


changing (EC) channel //^ //g, by using the extremely hig' 

fi~ —)■ e~ conversion experiments (COMET 


1 sensitivity of the ongoing 
40| and Mu2e [4^), very robust bounds have 


been set on the vector parameters j48|. To this end, we conclude that, if the Mu2e and 


COMET experiments wi 
parameters extracted in 
this channel. 


_not detect muon to electron conversion events, then the new ej^ 


48l | will lead to undetectable coherent rates at the SNS facility for 


For our present calculations we use the current bounds 48(] set by the sensitivity of the 
PSI experiment 8 ll and found countable number of events for the near detectors in the 
case of the corresponding //„ —)■ reaction. 


’he other exotic parameters, i.e. with 
a = e, and have been taken from Ref. [^. As discussed in Ref. [^, we do not take 


into account the contribution, since the corresponding limits are poorly constrained and 
eventually predict unacceptably high rates. 

Before closing, it is worth noting that, the present calculations indicate signihcant pos¬ 
sibility of detecting exotic neutrino-nucleus events through coherent scattering in the afore¬ 
mentioned experiments. Since neutrino-physics enters a precision era 8 |, a difference from 
the Standard Model predictions leads to an undoubtable evidence of non-standard neutrino- 
nucleus interactions (NSI). We recall that, in order to experimentally constrain simultane¬ 
ously all the exotic parameters at high precision, the detector material should consist of 
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T^ves. j.^res. 


FIG. 8. The expected non-standard neutrino scattering events over the recoil energy threshold 
at the COHERENT detector, filled with (left) 391 kg of liquid ^°Ne and (right) 456 kg of liquid 
^°Ar, both located at a distance of 20 m (^y . = 2.5 X 10^ i/s ^cm from the source. A perfectly 
efficient detector and negligible neutrino oscillation effects are assumed. 


maximally different ratio k 


(A + N)/(A + Z) l8|,l47|. 


Our Future plans include estimation of the incoherent channel which may provide a 
signihcant part of the total cross section, especially for energies higher than ^ 20 — 40 
MeV (depending on the nuclear target and the particle model predicting the exotic 
process). 
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FIG. 9. Same as Fig. [51 but for 100 kg of at 20 m = 2.5 x 10^ vs ^cm and 100 kg of 
liquid at 40 m . = 6.3 X 10 ® i/s ^cm from the source. 

V. SUMMARY AND CONCLUSIONS 

Initially, in this paper the evaluation of all required nuclear matrix elements, related to 
Standard Model and exotic neutral-current i/-nucleus processes is formulated, and realistic 
nuclear structure calculations of i/-nucleus cross sections for a set of interesting nuclear 
targets are performed. The hrst stage involves cross sections calculations for the dominant 
coherent channel in the range of incoming neutrino-energies 0 < < 150 MeV (it includes 

i/-energies of stopped pion-muon neutrino decay sources, supernova neutrinos, etc). 

Additionally, new results for the total number of events expected to be observed in one 
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ton of various z/-detector materials are provided and the potentiality of detecting supernova 
as well as laboratory neutrino-nucleus events is in detail explored. The calculations are 
concentrated on interesting nuclei, like ^'^Ne and ^°Ar, ^®Ge and ^^^Xe which are important 
detector materials for several rare event experiments, like the COHERENT at Oak Ridge 
National Laboratory, and also experiments searching for dark matter events as the CERDA, 
SuperCDMS, XENON 100, CLEAN, etc. By comparing our results with those of other 
methods, we see that the nuclear physics aspects (reflecting the accuracy of the required 
z/-nucleus cross sections), appreciably affect the coherent gs —)■ gs transition rate, a result 
especially useful for supernova z/-detection probes. 

In the present work, the QRPA method that considers realistic nuclear forces has been 
adopted in evaluating the nuclear form factors, for both categories of z/-nucleus processes, 
the conventional and the exotic ones. Also, a comparison with other simpler methods as (i) 
effective methods and (ii) the method of fractional occupation probabilities, which improves 
over the simple Shell-Model and gives higher reproducibility of the available experimental 
data, is presented and discussed. We conclude that among all the adopted methods the 
agreement is quite good, especially for light and medium nuclear isotopes. However, since 
coherent neutrino-nucleus scattering can probe the neutron nuclear form factors, methods 
like the BCS provide more reliable results. 

In view of the operation of extremely intensive neutrino fluxes (at the SNS, PSI, J- 
PARC, Fermilab, etc.), the sensitivity to search for new physics will be largely increased, and 
therefore, through coherent neutrino-nucleus scattering cross section measurements, several 
open questions (involving non-standard neutrino interactions, neutrino magnetic moment, 
sterile neutrino searches and others) may be answered. Towards this purpose, we have 
comprehensively studied the non-standard neutrino-nucleus processes and provided results 
for interesting nuclear detectors. Our predictions for the total number of events indicate that, 
within the current limits of the respective flavour violating parameters, the COHERENT 
experiment may come out with promising results on NSI. Moreover, this experiment in 
conjunction with the designed sensitive muon-to-electron conversion experiments (Mu2e, 
COMET), may offer signihcant contribution for understanding the fundamental nature of 
electroweak interactions in the leptonic sector and for constraining the parameters of beyond 
the SM Lagrangians. 
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